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Two unknown values among three electrochemical values, i.e. electrode area, diffusion coefficient, and concentration, 
are simultaneously obtained by nonlinear regression analysis of a single chronoamperometric faradaic current curve 
al a microdisk electrode. The approach is an analytical application of the semi-empirical equation presented by Shoup 
and Szabo for the chronoamperometric response at a disk electrode. To demonstrate the usefulness and accuracy 
of this approach, the chronoamperometric current at a platinum disk electrode of 50 urn radius in solutions of RufNHilc" . 
ferrocene, Fe(CN)j", and C M , wc:e analyzed. 



Introduction 

The chronoamperometric response at a disk electrode has 
been studied for determination of the diffusion coefficient 
of electroactive species 1_J . In most cases, at least two sepa- 
rate measurements and the approximate expressions for ch- 
ronoamperometric cun-enr were used. In this study, two un- 
known parameters among diffusion coefficient CD), concent- 
ration (C), and electrode radius (a), are simultaneously dete- 
rmined from a single ehronomperometric faradaic current 
at a disk electrode of small radius. 

In general, the physicochemical unknowns are obtained 
by the endeavor to fit the experimental data with the propos- 
ed model. Therefore, if the equation for the model is proven 
to be accurate, it can be applied to obtain unknown values. 
The diffusion controlled chronoamperometric response at a 
disk electrode was solved by Ksenzhek el al? and Aoki et 
atF. The former describes the current function as an integral 
equation of the 1st order Bessel function and the latter de- 
scribes that as a serial sum of Gamma function. The calcula- 
tion itself of these functions may be challenging. Therefore, 
both are hardly applicable to fit the experimental data. A 
simple semi-empiricat nonlinear equation which is accurate 
to 0.6% in arbitrary time domain was presented by Shoup 
and Szabo'. In the electrochemical experiment where 1% 
error is allowed, it may be regarded as a de facto solution. 



In principle, since the accurate equation is known, it can 
be applied to the experimental data for the analytical pur- 
pose. It can be easily analyzed by the nonlinear curve fitting 
or the nonlinear regression analysis using ubiquitous person- 
al computers in electrochemical laboratories, especially a co- 
mputer with a data acquisition board. In the potential step 
experiment, the charging current is not negligible at the be- 
ginning time, but quickly (exponentially) decreases with time', 
it. i, « exp {-l/(.R,C u )) where R, ie the solution resistance 
and C is the double layer capacitance. Therefore, the time 
domain for dominated faradaic current can be adjusted to 
fit the Shoup and Szabo equation which describes the fara- 
daic current at a disk electrode. The time domain can be 
extended to very short time for the high concentration of 
electroactive species and electrolytes. For the arbitrary con- 
centrated solutions, the chronoamperometric faradaic current 
should be obtained by the total current subtracted by non- 
faradaic current which can be approximated from the current 
in the same electrolyte solution without the electroactive 
species and in the sarne experimental conditions. In this 
work, the latter method is employed. 

To demonstrate the usefulness and accuracy of the nonli- 
near regression approach, it is applied to the chronoampero- 
metric faradaic current at a platinum disk electrode of 50 
urn radius in various solutions. In result, we report values 
of diffusion coefficients of ferricyanide in aqueous solution, 
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rutheniurndlD hexaammine a( pH 7.0 in aqueous solution, 
ferrocene in acelonitrile, and Cm in dichloromeihane. The 
major goal in this study is the demonstration of the nonlinear 
regression approach by which two unknown values among 
three electrochemical values, ie. electrode area, diffusion 
coefficient, and concentration, can be simultaneously deter- 
mined when one of them is known. Note that the » value 
of electrode reaction, O + ne -»R, is assumed to be known 
throughout this work. 

Theory 

Considering the semi-infinite spherical diffusion for the 
electrode reaction, O+ne -»R, the principle of a new meth- 
od can be easily understood because the equation is simple 
and linear. The chronoamerometric current* 1 at a sphere elec- 
trode is as follows. 

»(/) = nFADC/r. + nFADCVfnDty 1 

= nf4nr.DC* 4 tiFi^r.WCV " « (1) 

, where A is the area of electrode, n is the number of elec- 
trons transferred per molecule, F is the Faraday constant, 
r. is the radius of spherical electrode, D is the diffusion 
coefficenl, C* is the bulk concentration of the electroactive 
species, and t is the elapsed time since the potential step 
is applied. As mentioned above, the n value is assumed to 
be known. Jn equation (]), »'(;) is a linear function of t~ v *. 
A plot of current versus / M has the following slope and inter- 
cept: 

slope=«F4n'Vi? w C* 
intercept -nFA-nrDC 

Thus, two of n. r„, D. and C. can he simultaneously deter- 
mined from 2 single chronoamperometric cuive, where (wo 
of them is known. H. Ikeuch el a/. 1 " used this approach to 
determine the diffusion coefficients of serveral metal ions 
in aqueous solutions with a hanging mercury drop electrode. 

The chronoamperometrk response at a disk embedded in 
art infinite insulating plane is accurately (0.67o error) describ- 
ed by a simple semi-empirical equation which is presented 
by Shoup and Szabo. 

f = 4nFDC'a[A, +4*- '"-Mjexpf-A*- '•*)] (2) 

, where z-m/a 1 and a is the radius of disk electrode. A h 
A 2 .Aj, and At are constant numeric values, A, =n/4, A:=n ,/2 / 
2. v4 3 =l-7t/4, and ^=( n w /2-4n-*«)/(l-ji/4). 

The current at a disk embedded in an insulator is similar 
to that at a sphere electrode because both approach to a 
nonzero steady-state current at infinite time (or / l,7 ->0) 
while the theoretical steady-state current at a plane electrode 
or a shielded disk electrode should be zero (Figure 1(A)). 
The direct comparison of cunenls al a sphere and disk of 
same radius is unfair because the surface areas are different, 
By this reason, a quasi-hemisphere electrode of arbitrary 
solid angle, ft is imaginatively designed. Note that the solid 
angle of full sphere is 4n and that of hemisphere is 2n. The 
conditions to determine the solid angle are that the surface 
area (A) and the steady-state current (»'„) of a quasi-hemis- 
phere electrode are same as those of a disk electrode, U 




Figure 1. Chronoamperometric response at three different elec- 
trodes. i/{4nFDC*o) versus t _w . (A) Current at a shielded disk 
electrode of radius a or at a plane eletrode of na 1 area, oottrell 
equation". (B) Current at a disk electrode of radius a, Shoup 
and Szabo eqution'. (C) Current nl a quasi-hemisphere electrode 
of radius (rc/4)a and 16/n solid angle". 

A = Br. , =na i (3) 

i„=nFDC*$r.=inFDC*a (4) 

From the above two conditions, 

r. = W4)a (5) 

0=16/n (6) 

In result, the currents at a disk electrode (Figure 1(B)) and 
a quasi-hemisphere electrode uf 16/n (or 1.6211n) solid angle 
(Figure 1(C)) are very similar throughout all time domain 
with 6% relative deviation". 

Since the chronoamperometric response at a disk electrode 
is similar to that of a quasi-hemisphere electrode, where 
r,=(n/4)o = 0.78543, and solid angle = 16/n = 1.621 bt, the 
same strategy can be applied to the analysis of a disk elec- 
trode. One minor difficulty comes from that the theoretical 
current at a disk electrode is nonlinear against t~ m . A semi- 
empirical equation presented by Shoup and Szabo" is simple 
but still nonlinear. Thus, the chronoamperometric current 
at a disk electrode can not be analyzed from approach used 
by H. Ikeuchi el al. v> . In this case, it can be analyzed by 
the nonlinear regression analysis. Theoretically, several pa- 
rameters can be obtained from the nonlinear regression ana- 
lysis. But, since the equation (2) closely resembles the cur- 
rent at a quasi-hemisphere electrode, practically two parame- 
ters can be obtained. For the nonlinear regression analysis, 
the algorithm, which is well described in the textbooks" 13 , 
is employed. The computer coding of this algorithm Is con- 
structed by C-language and compiled by Tiirbo-C 3.0. Since 
the procedure is well described in refences 12 and 13. the 
detailed description is omitted. In brief, the probleuvis to 
find the best values of D, C* and a, to fit the experimental 
data curve (current versus time) in the appropriate time do- 
main. Starting with the initially guessed values of D, C*, 
and a, the better guessed values to minimize the chi-square 
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merit function" are found. The procedure is repealed until 
cWi square merit (unction stop* decreasing. 

Experimental 

AH chemicals were used as received from commercial sou- 
rces except for C w which was obtained by the purification 
of soot produced in plasma discharge of graphite rod" 
The supporting electrolyte of ferrocene (Aldrich Chemical Co.. 
Inc.) was 0.1 M tetraethytanunonium perchlorate (BAS, Inc.) 
in acetonitnle (Merck). For potassium ferricyanide (Sigma 
Chemical Co.), the supporting electrolyte was 1.0 M KCI (Aid- 
rich Chemical Co.. Inc.). The supporting electrolyte for ru- 
thenium hexaamine trichloride (Aldrich Chemical Co., Inc.) 
was 0.1 M Na?HPO< (Junsei Chemical Co., Ltd.) adjusted 
to pH 7.0 with phosphoric acid. For Co. the supporting elec- 
trolyte was 0.1 M tetrabutylammonium tetrafluoroborate (Al- 
drich Chemical Co., Inc.) in dichloromethane (Junsei Chemi- 
cal Co., Ltd.). All solutions were deoxygenaled by argon bub- 
bling prior to measurements. Aqueous solutions were prepa- 
red with triply distilled water. 

A conventional three electrode system with one compart- 
ment cell was used for all measurements. A platinum wire 
(0.5 mm diameter and ca. 5 cm length) was used as a counter 
electrode. A platinum disk electrode of 50 urn radius (BAS, 
Inc.) was used as a working electrode. As a reference elec- 
trode, a saturated calomel reference electrode(SCE) was 
used for aqueous systems and a silver quasi-reference elect- 
rode (Ag-QRE) was used for nonaqueous systems. In this 
work, electrode was polished with 0.05 urn alumina solution 
or sequentially 6 um, 3 urn, and 1 um diamond paste, and 
used in 5 to 10 chronoamperometric experiments. 

For the cyclic vutammelry and chronoamperometry, a home- 
made polentiostat, which is interfaced to a personal computer 
(NCS-386S System. Newcom System Co., Korea) via a data 
acquisition card (12 bit A/D converter card, Model PCL-812, 
Advantech Co., Ltd., Taiwan), was used. The control program 
was constructed by C-language. 

Results and Discussion 

As mentioned above, the current at a disk (Figure 1(B)) 
is similar to a quasi-hemishpere electrode (Figure 1(C)). Both 
curves approach to nonzero intercepts at infinite Lime, while 
the current at a shielded disk electrode or planar electrode 
should approach to zero value. The nonlinearity of the cur- 
rent at a disk electrode can be exaggerated in the plot of 
dWitnFD&aWdfr-™) versus logd" 1 *) as seen in Figure 2. 

The chronoamperometric faradaic current was obtained by 
the total current subtracted by non-faradaic current which 
can be approximated from the current in the same electro- 
lyte solution without the elctroactive species and in the same 
experimental conditions. For the time domain to be analyzed, 
in general, the following two factors should be considered. 
One is the instrumental error and the other is the size of 
electrode. First of all, because of the data acquisition system 
of 5 KHz in this work, the time region of less than 1 msec 
is not possible. Also, the RC time constant in the current- 
voltage converter (100 nA/Volt) forces to discard data at 
short rime. In this work, the fixed feedback system of the 
fixed RC time constant is measured as 9<± 1) msec using 
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Figure 2. A plot of d(i/(inFDC m a))Jd[z- ,!l ) vs. log (r "*) for the 
chronoamperometric response. (A) Current at a shielded disk 
elertrode of radius a. (B) Current at a disk electrode of radius 
a. (C) Current at a qua.si-hemispherc electrode of radius (n/4)a 
and 16/71 solid angle. 

a dummy cell. Therefore, the current data after the 5 times 
of this is in less than 1% error". In result, the shortest time 
with this instrument is 50 msec On the other hand, if the 
current data are acquired only at the longer time region, 
then the variations of current and l xn is so small to be 
alike a point data in the plot of i versus / m . In numerous 
trials, it is found that the optimum time domain for the non- 
linear regression analysis is in r" ,/I — 03 to 2.5. when the 
current curve is distinctively nonlinear as seen in Figure 
2. But, the data before 50 msec has to be rejected to avoid 
the instrumental error as mentioned above. For example, 
when the radius of disk is 50 um and the diffusion coefficient 
is 0.5X10" 5 cmVsec as a typical value, the real time domain 
is in the range of 0.2 to 14 sec But, in the case of 5 um 
radius, the real time domain is in the range of 2 to 140 
msec The current voltage converter of 10 times higher sen- 
sitivity with decreased RC time constant is mandatory. The- 
refore, the platinum disk electrode of 50 urn radius was em- 
ployed in this work. The larger disk electrode of 1 mm dia- 
meter can be employed. In this case, the vibrational problem 
has to be considered, because the time domain has to be 
extended to 1400 sec without disturbing the larger diffusion 
length {2Dt)" 2 =L2 mm. 

The half-wave potential is determined by cyclic vollamme- 
tric curves as listed in Table 1. For the diffusion controlled 
chronoamperometry, the initial potential is selected as far 
negative from this half-wave potential for oxidation (or far 
positive for reduction) and the step potential is selected as 
far positive for oxidation (or far negative for reduction). 
These potential values are listed in Table L In this way, 
the step potentials can be unambiguously selected for the 
organic solvent system where the unreliable Ag-QRE is used. 

The typical curve fitting results of the nonlinear regression 
analysis are shown in Figure 3. All the experimental results 
are excellently fitted to the theoretical curve. The numerical 
results of nonlinear regression for Ru(NH 3 )J', ferrocene, and 
Fe(CN)J' are summarized in Table 2. First of all, the radius 
values are not much deviated from the specification value 
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Table 1. Hall-wave Potentials for Four Electroaclive Solutions 
and the Applied Potentials in Chronoamperometrk Measureme- 
nts 



Solutions' 


ftatmV) 


£.im\'t 


E.(nW 


Re(erenct 
electrode 


(1) RuCNKOj' 


-225' 


0 


-400 


SCE 


(2) Ferrocene 


55CK 


300 


800 


Ag-QRE 


(3) Fe(CN)3 


230 


500 


0 


SCE 


(4) Cm 


340 


-500 


-600 


Ag-QRE 



"(1) RufNHJg': 1.00 mM Ru(NJ-UA in the aqueous solution 
of 0.1 M NajHPOj. adjusted to pH=7.0 with phosphoric acid. 
(2) Ferrocene: 1.00 mM ferrocene in acentonitrile solution of 
0.1 M letraethylammanium perchlorate. (3) Fe(CN)J ': 4.00 mM 
KJetCNlc in the aqueous solution of 1.0 M KCL (4) C H ,: Cm 
in dichloromethane solution of 0.1 M letrabutylammnniurn tetraf- 
luoroborate. »£ is initial potential before potential step lo 
Scan rate = 305 mV/s 'Scan rale = 610 mV/s. 



Table 2. The Diffusion Coefficients CD) and Electrode Radii (<z) 
Determined Simultaneously from Nonlinear Regression Analysis. 
Every Reported Result is the Average Value CD or a) of the 
Nonlinear Regression Analysis Results of at Least 5 Repetitive 
Chronnamperometric Curves with a Fresh Working Electrode 
Which Was Polished with Alumina Solution or Diamond Paste. 
Units : 10' s cmVsec for D and 1 urn for a 



Solution* 
Experiment 


RuCNIUi!'' 


Ferrocene" 


Fe(CN)i-« 


D a 


D a 


D a 


1 


0.566 505 


2.46 51,5 


0.780 50.4 


2 


0.580 49.9 


2.49 51.7 


0.764 50.7 


3 


0.586 50.4 


2.45 515 


0.782 49.8 


Average 


(1.577 503 


2.47 51.6 


0.775 503 



"See Table 1. * Temperature- 28+ 0SC. Polished with 0.05 urn 
alumina solution. ' Temperature = 25+ 0.5t. Polished with 0.05 
um alumina solution. ''TemperaLure-27± OSt . Polished seque- 
ntially with fi um, 3 tim, 1 um diamond paste. 



(A) 



Figure 3. Comparison of experimental current (•> and therore- 
tical current (-). Ca) Ru(NH : ,) r r . (b) ferrocene, (c) Ke(CN)J". and 
(d) C». The sampling rate was 40 msec per point except for 
ferrocene where it was 20 msec. See Table 1 for the experimen- 
tal conditions, i.e. solvent and electrolyte. 



-1/2. 



of BAS, Inc., 50 um. The spread of data are less than ±.0.9 
um for one electrochemical system, Le. 49.9 to 505 um in 
RufNHJiT system, 51,5 to 51.7 um in ferrocene system, and 
49^to 50.7 um in Fe(CNfi' system. The radius measure- 



lationship" was considered, it the diffusion coefficient in- 
creases by 1 to 2% with increasing temperature. The diffusion 
coefficient value (0=0577X10 5 cmVsec at 28t ) of RuCNHjJT 
in Table 2 shows good agreement with that reported by Wight- 
man et al? CD=0-54SX10-' cmVsec at 25t ), assuming D 
value increases by 1.7* with temperature. The diffusion coe- 
fficient value {£>= 2.47X10 s cmVsec at 25'C) of ferrocene 
in Table 2 also shows good agreement with the literature 
value CD=2.4X10- & cmVsec at 22t>", assuming D value 
increases by 1.0% with temperature. In the case of ferricya- 
nide, the literature values of ferricyanide are scattered in 
broad range". The diffusion coefficient value of Fe(CN>2" 
in Table 2 is in good agreement with those reported by von 
Stackelberg CD=0.763X10- 5 cmVsec at 25t )'• and Adams 
(D~ 0.74 X 10- 5 cmVsec at 25t )», assuming D value inrreasfw 
by 0.8 or 2.3% with temperature, respectively. The former 18 
reported the resulting same D values employing two differ- 
ent electrodes. One is the carefully constructed planar plati- 
num electrode system which surely excludes the edge effect 
and the other was the dropping mercury electrode system 
(polarography). The latter* used the flow Injection analysis 
technique which was non-electrochemical method. The D val- 
ues of various methods, including the method in this -work, 
are all in good agreements. But, the value CD =0.726X10" 5 
cmVsec at 25t ) obtained by thin-layer cell 1 * seems measu- 
red low. 

In the case of C» solution in dichloromethane, three para- 
meter D, C*. and a, are assumed to be all uncertain. To 
solve this problem, two step method is used. At first, the 
radius of electrode is determined from the nonlinear regres- 
sion analysis of ferrocene solution system. The same elec- 
trode is immediately transferred to Ceo solution system to 
measure the chronoamperometric response. In this process, 
since the effective radius of electrode may charge as a result 
of polishing, the electrode is not polished at all. Instead it 
-was just rinsed with acetonitrile and dichloromethane and 
dried in the air. By this procedure, the results of nonlinear 
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Table 3. The Simultaneous Determined Diffusion Coefficient 
(D) and Concentration (C*) of C» from Nonlinear Regression 
Analysis. Every Reported Result is the Value CD, C* or a) Calcu- 
lated from the Nonlinear Regression Analysis of a Single Chro- 
noamperometric Faradaic Curve. A Working Electrode Was Po- 
lished with 0.05 urn Alumina Solution. Units : 10 " 5 cmVsec for 
D, 1 pm for a and 1 mM for C 

Solution" Ferrocene* CV CV C«' 

Experiment D a D C J) C* D C 

1 2.52 51.4 1.02 0219 0.98 0.232 1.07 0.207 

2 2.51 51.3 1.05 0.215 L01 0.228 1.09 0.203 
T 951 513 107 (V9.17 1 03 0.230 1 U O'W 

4 2.51 51.3 1.03 0.226 0.99 0.239 1.07 0.213 

5 2.51 51.2 1.12 0.214 1.23 0.201 1.16 0.202 
Average 2S\ 51.3 1.06 0.21B 1.05 0.226 1.10 0.206 

'See Table 1. 4 Temperature=-27± 0,5'C . 'Temperature = 25+ 1 
t. Assuming a =51.3 pm. < Temperature=25i it . Assuming 
pm. 'Temperature=25± It . Assuming a=52.3 jun. 

urn in 95% confidence level). But, it has to be reminded 
that the error associated with the determination of electrode 
radius between two different systems is 2% as mentioned 
above. The radius of electrode is a = 513± 1.0 pm. Using 
the minimum and maximum value of a Z>=(1.05± 0.12)X10 _S 
cmVsec and D=(1.10±0.05)X10- 5 cmVsec in 95% confide- 
nce level. Therefore, in 95% confidence level, D value of 
C« is (1.05+0.12)X10- 5 em'/sec. Diffusion coefficient of Cm 
in this work is about 2 times higher than that determined 
by Kadish el al. 21 using the rotating disk technique at 22± I 
t . However, considering the accuracy of this approach, the 
diffusion coefficient reported by Kadish el al. seems to be 
measured lower. 

Conclusion 

Two of the three parameters D. C*. and a. are simulta- 
neously determined from a single chronoamperometric fara- 
daic current curve at a microdisk electrode. It is demonstrat- 
ed that this approach can be consistently applied to any elec- 
trochemical systems with less than 2% error. This approach 
can be helpful in numerous analytical situations, especially 
when the standard curve cannot be obtained. Even when 
all the electrochemical parameters are unknown, a couple 
of chronoamperometric experiments lead to the determinate 
results of D. C* and a. Most importantly, the diffusion coeffi- 
cient can be precisely obtained without separate determina- 
tion of the effective surface area which has to be carefully 
obtained by the system of precisely known D and C* in 
the other electroanalytical method. In result, we reported 
the diffusion coefficients nf four electroactive species. Among 
them, only the diffusion coefficient of Cm in dichlorornethane 
shows the huge difference with the literature value, while 



the other three are consistently in good agreements with 
the literature values. Therefore, we conclude that the diffu- 
sion coefficient of Cm obtained in this work is more reliable. 
The diffusion coefficient of C w in dichloromelhane is (1.05± 
0.12)X10" S cmVsec in 95% confidence level. 
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